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Introduction 

Freak waves are very large water waves whose heights exceed the significant wave 
height of a measured wave train by a factor of more than 2.2. However, this in 
itself is not a well established definition of a freak wave. The mechanism of freak 
wave generation in reality as well as modeling it in a wave basin has become an 
issue of great importance. 

Recently one is aware of the generation of freak wave through the Benjamin-Feir 
type of instability or self focussing. Consequently the Non-Lincar-Schrodinger 
(NLS) equation forms a good basis for understanding the formation of freak 
waves. However, the complex generation of a freak wave in nature within a sea 
condition is still not well understood, when the non-linearity of the carrier wave 
is not small. In our study we will focus on the Soliton on Finite Background, 
an exact solution of the NLS equation, as a generating mechanism for extreme 
waves. 

Apart from a numerical investigation into the evolution of a soliton on a fi- 
nite background also extensive detailed model tests have been performed for 
validation purposes in the hydrodynamic laboratories of the Maritime Research 
Institute Netherlands (marin). Furthermore, a numerical wave tank [12] is used 
to model the complete non-linear non-breaking wave evolution in the basin. 

Properties of the Soliton on Finite Background 

The NLS equation is chosen as a mathematical model for the non-linear evolution 
of the envelope of surface wave packets. For spatial evolution problems, it is given 
in non-dimensional form and in a frame of reference moving with the group 
velocity by 

d £ V> + i/3<9rV> + n|V'lV = 0, (1) 

where £ and r are the corresponding spatial and temporal variables, respectively; 
(3 and 7 are the dispersion and non-linearity coefficients. This equation has many 
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families of exact solutions. One family of exact solutions is known as the Soliton 
on Finite Background (SFB) and this is a good candidate for describing extreme 
waves. This exact solution has been found by Akhmediev, Eleonskii & Kulagin 
[3J, see also [2] and pQ. 

This SFB solution describes the dynamic evolution of an unstable modu- 
lation process, with dimcnsionless modulation frequency v. In the context of 
water waves, for infinitesimal modulational perturbations to a finite-amplitude 
plane wave, this process is known as Benjamin-Feir (BF) instability [5]. How- 
ever, non-linearity will limit this exponential growth and the SFB is one (of 
many other) non-linear extension of the BF instability for larger amplitudes of 
the modulation. Extensive research on the NLS equation and the SFB solution, 
to obtain a better understanding of deterministic extreme-wave phenomena has 
been conducted in the past few years (see e.g. [TO], [H], [I] and [OJ). 

An explicit expression for the SFB is given as the following complex-valued 
function 

i ft ~ \ J v 2 cosh(<7Q - i [<r/(7rg)] sinh(gg) \ 
ip{£,T;v,r )=A{€)-i L 1>, (2) 

[ cosh(<7£) — i I — \v 2 cos(^t) J 

where A(£) — r e~ nr °£ is the plane-wave or the continuous wave solution of the 
NLS equation, a = ^r^v — v 2 is the growth rate corresponds to the Benjamin- 
Feir instability, v = i>r Q ^J^ is the modulation frequency, and v, < v < y/2 is 

the normalized modulation frequency. This SFB reaches its maxima at (£, r) = 
(0, ^t 5 -), 7i G Z. It has a soliton-like form with a finite background in the spatial 
^-direction. The SFB is periodic along in the temporal r-direction, with period 
=^-. For |^| — > oo, the SFB turns into the continuous wave solution A(£). It 
possesses two essential parameters: r n and v. 
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Fig. 1. Density plots of a physical wave packet profile according to an SFB envelope 
for ii\ = 1, showing the wave dislocation phenomenon. 
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The first-order part of the corresponding physical wave packet profile rj(x, t) 
for a given complex-valued function ?/;(£, r) is expressed as follows 

r)(x, t) = T y^-^ot) + C C t ( 3 ) 

where c.c. means the complex conjugate of the preceding term, the wave number 
fco and frequency uiq satisfy the linear dispersion relation ui — f2(k) = V^'tanh/c. 
The variables (x, t) in the non-moving frame of reference are related to (£, r) in 
the moving frame of reference by the transformation £ = x and t = t — x/f2'(ka). 
The modulus of "0 represents the wave group envelope, enclosing the wave packet 
profile r](x,t). The dimensional laboratory quantities are related to the non- 
dimensional quantities by the following Froude scaling, using the gravitional 
acceleration g and the depth of the basin h: x\ a ^ — x ■ h, t\ a h = t ■ v/f , &i a b = k/h, 
^iab = w ■ yf, and 7/ lab =r)-h. 

In principle, the wave profile including the higher-order terms represents a 
good approximation to the situation in real life. To accommodate this fact, we 
will include higher-order terms up to second order. We apply an perturbation- 
series expansion (Stokes' expansion) to the physical wave-packet profile r](x,t) 
and the multiple-scale approach using the variables £ and r, where £ = e 2 x, 
t = e(t — x/i?'(fco)), and e is a small positive non-linearity and modulation 
parameter. The corresponding physical wave packet profile, consisting of the 
superposition of the first-order harmonic term of 0(e) and a second-order non- 
harmonic long wave as well as a second-order double-frequency harmonic term 
of 0(e 2 ), is given by 

n(x, t) = e U (10) (£, tY^*-^ + c. 

+ e 2 { V (20) & r) + [V (22) (£, T )e 2i ^ x -^ + c.c] } . (4) 

We find from the multiple-scales perturbation-series approach that ip^H^ T ) = 
?H£,t) satisfies the spatial NLS equation and 

M(t rl L_ ^Wfco) - ft(*0) i .(10W rl |2 f5l 

^ (22) fer) = fco '" ta "3?° [^^fcr)] 2 . (6) 
2 tanh Ko 

A similar derivation for the temporal NLS equation resulting from the KdV 
equation can be found in [B], By including this second-order term, the wave signal 
r)(x,t) experiences the well-known Stokes' effect: the crests become steeper and 
the troughs becomes shallower [§]. 

The coefficients j3 and 7 of the spatial NLS equation are given, in nondimen- 
sional form, as: 

P = - 1 (7) 
' 2[r2'(fc )] 3 ' y> 

7= W ' (8) 



4 



where 



, 2 9 tanh 4 k - 10 tanh 2 k + 9 



7i = k 2 Q f2(ko) 77-^47 — > (9) 

4 tanii fco 

1 2 l-tanh 2 fc 

A= 2 fc " r2(fc ) ' (10) 

1 4fcofl'(fc ) - g(fcp) 
" c Q{k ) [Q'{Q)f - [Q'{k )Y and 

^ = «C^o)-^y d2) 
These can be used to compute the SFB solution r) from Equation |2]). 



Phase singularity and wave dislocation 




Fig. 2. The evolution of the SFB for v\ = \ from a modulated continuous wave signal 
into the extreme position. From top to bottom, the signals are taken at x = —200, 
x — —100, x — —50, x — —30, and x — 0. 
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By writing the complex-valued function ip in a polar (or phase-amplitude) 
representation, it is found that for modulation frequencie v in the range < v < 

a phase singularity phenomenon occurs. It happens when the real-valued 

amplitude vanishes and therefore there is no way of ascribing a value to the 
real-valued phase when it occurs. The local wave number k = ko + d x 8 and local 
frequency u) = ujQ — dtO, with r) = arg (^(C, t )j become unbounded when this 
happens. The corresponding physical wave packet profile rj(x, t) confirms this by 
showing a wave dislocation phenomenon. When the real- valued amplitude \tp\ 

vanishes at that specific position and time, waves merge or split. For < 

v < \/2, the real-valued amplitude is always definite positive, and thus there 
is no wave dislocation. Furthermore, in one modulation period, there is a pair 
of wave dislocations. Before or after this dislocation, the real-valued amplitude 
reaches its maximum value. Figure [l] shows the density plot of a physical wave 
packet profile 77(2, t). The wave dislocation is also visible in this figure. Figure[2] 
shows the evolution of the SFB from a modulated wave signal until it reaches the 
extreme position. We can see also in this figure that the amplitude \ip\ vanishes 
at some moments for the extremal position x — 0, causing phase singularity. 

The phase singularity is a well known phenomenon in physical optics. In 
the context of water waves, similar observations can be made, and also wave 
dislocations occur. Trulsen [13] calls it as crest pairing and crest splitting and 
he explains this phenomenon as a consequence of linear dispersion. 



Maximum temporal amplitude 

The maximum temporal amplitude (MTA) is a useful concept to understand 
long-time behavior of wave elevation. For wave propagation in the laboratory, 
it also gives a direct view of the consequences of an initial wave signal on the 
corresponding extreme-wave signal. It is defined as 

[i(x) = max t rj(x, t), (13) 

where r](x,t) is the surface elevation as a function of space x and time t. It 
describes the largest wave elevation that can appear at a certain position. For 
laboratory wave generation, it describes the boundary between the wet and dry 
parts of the wall of the basin after a long time of wave evolution. 

Figure|3]shows the MTA plot of the SFB in the laboratory coordinates. In this 
example, the mean water depth is 3.55 m and the wavelength is approximately 
6.2 m. The wave signal is generated at the left side, for example at xub = —350 
m, and it propagates to the right and reaches its extremal condition at xuh = 0. 
A slightly modulated wave train increases in amplitude as the SFB waves travels 
in the positive x-direction. Furthermore, in this example a SFB wave signal with 
initial amplitude around 0.19 m can reach an extreme amplitude of 0.45 m, an 
amplification factor of around 2.4. After reaching its maximum amplitude, the 
MTA decreases monotonically and returns to its initial value. 
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Fig. 3. The MTA plot of the SFB, the corresponding wave profiles at t = 0, and its 
envelope. 



Experimental Result 

For the validation of the proposed method we performed experiments in one of 
the wave basins of marin. The basin dimensions amounted to L x B x D as 
200 m x 4.0 m x 3.55 m. In the basin an array of wave probes were mounted as 
indicated in the set-up in figure Q. The predefined wave board control signal 
was put onto the hydraulic wave generator. The stroke of the wave flap was 
measured. Main characteristics of the model test experiments: 



200m 

















1 OOm 











Fig. 4. The set-up of the wave probe array in the wave basin 



Carrier wave period is 1.685 sec. maximum wave height to be achieved (MTA) 
varies from 0.213m to 0.2485m 
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As a explanation the results for the tests with an MTA of 0.2485 will be shown, 
see figure ([5]) to figure ([7]). 
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Fig. 5. Comparison Non-Linear Wave model HUBRIS with results from experiments 
and sNLS 

Figure [8] shows the SFB signal based on the experiment at distance 150 m 
from the wave maker, where it is expected that the signal to be extreme. That 
figure also shows the phase singularity phenomenon when the local frequency 
becomes unbounded when the real-valued amplitude vanishes or almost vanishes. 
The experiment result shows asymmetric form of the extreme signal while the 
theoretical result of the SFB preserve the symmetry of the signal. It is suspected 
that if the modified NLS equation of Dysthe [7] is used as the governing equation 
for the wave signal evolution, then there are good comparisons with experimental 
measurements. The good comparisons are observed for the case of bi-chromatic 
waves, where the modified NLS equation predicts both the evolution of individual 
wave crests and the modulation of the envelope over longer fetch [T3] . 

References 

1. Akhmediev, N.N. & Ankiewicz, A. "Solitons, Nonlinear Pulses and Beams", 
Optical and Quantum Electronic Ser., Vol. 5, Chapman & Hall, l s ^ edition, (1997). 




150 160 170 180 190 2 210 2 2 23 240 250 

t Tsl 
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Fig. 7. Comparison Non-Linear Wave model HUBRIS with results from experiments 
and sNLS 
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Extreme signal, x = 1 50 m, M = 28.73 cm, with predicted M = 21 .30 cm, v = 1 
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Local frequency, x = 1 50 m, M = 21 .30 cm, v = 1 
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Fig. 8. The SFB signal plot based on the experiment at 150 m (top) and the corre- 
sponding local frequency plot (bottom). 



